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Abstract 

We revisit the one-loop calculation of the annihilation of a pair of the lightest neutralinos 
into a pair of photons, a pair of gluons and also a Z photon final state. For the latter we have 
identified a new contribution that may not always be negligible. For all three processes we 
have conducted a tuned comparison with previous calculations for some characteristic scenarios. 
The approach to the very heavy higgsino and wino is studied and we argue how the full one- 
loop calculation should be matched into a more complete treatment that was presented recently 
for these extreme regimes. We also give a short description of the code that we exploited 
for the automatic calculation of one-loop cross sections in the MSSM that could apply both 
for observables at the colliders and for astrophysics or relic density calculations. In particular 
the automatic treatment of zero Gram determinants which appear in the latter applications is 
outlined. We also point out how generalised non-linear gauge fixing constraints can be exploited. 



LAPTH-1106/05 



tUMR 5108 du CNRS, associee a I'Universite de Savoie. 



1 Introduction 



We now have overwhelming evidence that ordinary matter accounts for a minute portion 
of what constitutes the Universe at large. Most impressive is the confirmation from the 
very recent WMAP datapp. Very interestingly, many extensions of the standard model 
whose primary aim was related to the Higgs sector and the mechanism of symmetry 
breaking do provide a good dark matter candidate. Very soon, with the energy frontier 
that will open up at the LHC, intense searches for this new physics with its associated 
dark matter candidate will be pursued in earnest. 

Meanwhile, many astroparticle experiments are going on, and will be improved by the 
time the LHC runs, to detect dark matter particles. The problem, either for direct or 
indirect detection of dark matter outside the colliders, is that we do not control many 
astrophysical parameters. For indirect detection which is the result of the annihilation of 
a dark matter pair in, say, the galactic halo of our galaxy, the photon signal is cleaner 
than that of the charged positron and antiproton that are considered as sources of exotic 
cosmic rays. The photon will point back to the source while the antiproton flux suffers 
from uncertainties due to the propagation. Of course in both cases one still has to rely 
on a modelling of the dark matter profile since one needs to know the number density of 
the annihilating dark matter particles. A very distinctive signal though would be that 
of a "direct" annihilation into a monochromatic photon. In this case the spectrum will 
reveal a peak at an energy corresponding to the mass of the annihilating particles since 
the latter move at essentially zero relative velocity v. In the galactic halo vjc ^ 10^^. 
Therefore, provided one has a detector with good energy resolution, the flux from the 
"direct" annihilation will clearly stand out above the (astrophysical) background or the 
diffuse contribution. The latter is due essentially to annihilation into quarks and W 
which subsequently fragment and radiate/decay into photons. This contribution has a 
continuous featureless energy distribution which is only cut-off at a maximum energy 
corresponding to the mass of the dark matter particle. There are, and there will be, 
many powerful detectors to search for such photon signals, covering a wide range in 
energy from MeV to TeV. These are either ground based, like the atmospheric Cerenkov 
telescopes, ACT,(Cangaroo^, HESS[3j, MAGIC[4j, VERITAS [5^,..) or space borne telescopes 
(EGRET[7J, AMS[6:, the upcoming GLAST[S],..). Many see in some of the present data an 
excess that might be a sign for New Physics and dark matter annihilation but we should 
probably be cautious and await confirmation from other more precise detectors covering 
the same energy range. One should also improve on the theoretical predictions and a 
better understanding of the background and the astrophysical component that enter the 
calculation of the photon yield. 

Our aim in this paper is to revisit the calculation of the "direct" self-annihilation into 
77, Z7 and gg of the lightest supersymmetric particle (LSP). This is a neutralino that 
we will denote as Xi- There have been a few attempts of calculating the one-loop induced 
XiXi ~^ 11 before two complete calculations jSJ settled the issue. These calculations have 
been made in the limit f = as is appropriate for dark matter annihilation in the halo. 
A very recent calculation[10j has also been made for this mode. Their results for u = 
agree in their most important features (higgsino limit, for example) with those of Refs jH], 
but as far as we are aware no systematic comparison has been performed. Much more 
important however is that there is, at the moment, only one calculation of XiXi ~^ ^1 [TT] 
foerformed at f = 0) desDite the fact that new features aDDear in this comoutation. These 



features, as we will see, can not be a mere generalisation of the 27 final state. We will 
in this paper calculate both XiXi ~^ 77 ^"^^ XiXi ^7 fo^' velocity and make a 
tuned comparison with the existing codes for v = 0, DarkSUSYjl2j for XiXi ~^ 77 
XiXi ~^ and PLATONdml for XiXi ~^ 11 ■ Iii XiXi ~^ ^7 we have identified a new 
contribution not taken into account in^I]. We will also show some results for v = 0.5 for 
both processes. 

As a by-product we will also compute the self-annihilation into gluons: XiXi ~^ 99 CH]- 
This can be derived from XiXi ~^ 77 by only keeping the coloured particles and dealing 
properly with the colour structure. This process could contribute to, for example, the 
antiproton signal. We will see that our results for XiXi ~^ 99 completely agree with those 
of DarkSUSY[T2] and PLATONdmg[10, for = and with PLATONgrelfTHj for v = 0.5. 

As is known pilTT|. the largest contributions to XiXi ~^ 77 XiXi ~^ ^7 ; especially for 
large neutralino masses, occur when the neutralino is a wino or a higgsino. As first pointed 
out in E] the cross section times the relative velocity, av, for both modes, tends to 
an asymptotic constant value that scales as for v = and large LSP mass. This 

result which breaks unitarity is due to the one-loop treatment of a "threshold" singularity 
that is nonetheless regulated by M^. It has very recently been, admirably, shown [Ti] 
how to include the higher order corrections through a non-relativistic non-perturbative 
approach. The latter reveals the formation of bound states with zero binding energy that 
show up as sharp resonances that dramatically enhance the cross section for particular 
masses. We have therefore thought it worthwhile to study the one-loop derivation in these 
scenarios and see how one can match the non-perturbative regime. The reason we do this 
and the main reason we carry the calculation of XiXi ~^ 11 ^^id XiXi ~* ^1 is that one 
needs a reliable code for the photon flux from self-annihilating neutralino LSP's. These 
cross sections have been missing from micrOMEGAs[15 that we have been developing for 
a very accurate derivation of the relic density in supersymmetry and which is currently 
being adapted to direct and indirect detection. The present paper only deals with the 
cross section calculation, leaving aside the astrophysical issues to the implementation 
and exploitation within micrOMEGAs. See however Ref. |^ for a preliminary use of 
micrOMEGAs to indirect detection using some of the results of this paper. 
The results presented in this paper constitute some of the first applications of a code 
for the automatic computation of one-loop processes in supersymmetry relevant both 
for the colliders and astrophysics, such as the problem at hand. Most crucial for the 
latter is a careful treatment of the loop integrals since for these applications the use of 
general libraries is not appropriate leading to division by zero because of the appearance 
of vanishing Gram determinants in the reduction of the tensor integrals. We will show 
how to easily circumvent this problem. 

2 Set-up of the automatic calculation 

Even in the standard model, one-loop calculations of 2 — > 2 processes involve hundreds of 
diagrams and a hand calculation is practically impracticable. Efficient automatic codes 
for any generic 2^2 processes, that have now been exploited for many 2 — > 3^71 ^] 
and even some 2 4^21120] processes, are almost unavoidable for such calculations. For 
the electroweak theory these are the GRACE- loop [21] code and the package FormCalcj22j 



based on FeynArts|23j and LoopToolsj 

With its much larger particle content, far greater number of parameters and more complex 
structure, the need for an automatic code at one-loop for the minimal supersymmetric 
standard model is even more of a must. A few parts that are needed for such a code 
have been developed based on the package FeynArtsusy[25^ but, as far as we know, no 
complete code exists or is, at least publicly, available. Grace-susy[26 is now also being 
developed at one- loop and many results exist |27j. One of the main difficulties that has 
to be tackled is the implementation of the model file, since this requires that one enters 
the thousands of vertices that define the Feynman rules. On the theory side a proper 
renormalisation scheme needs to be set up, which then means extending many of these 
rules to include counterterms. When this is done one can just use, or hope to use, the 
machinery developed for the SM, in particular the symbolic manipulation part and most 
importantly the loop integral routines and tensor reduction algorithms. 
The calculations that we are reporting here are based on a new automatic tool that uses 
and adapts modules, many of which, but not all, are part of other codes. We will report on 
this approach elsewhere. Here we will be brief. In this application we combine LANHEP 
(originally part of the package C0MPHEP[29 ) with the FormCalc package but with an 
extended and adapted LoopTools. LANHEP is a very powerful routine that automatically 
generates all the sets of Feynman rules of a given model, the latter being defined in a 
simple and compact format very similar to the canonical coordinate representation. Use 
of multiplets and the superpotential is built-in to minimize human error. The ghost 
Lagrangian is derived directly from the BRST transformations. The LANHEP module also 
allows to shift fields and parameters and thus generates counterterms most efficiently. 
Understandably the LANHEP output file must be in the format of the model file of the 
code it is interfaced with. In the case of FeynArts both the generic (Lorentz structure) 
and classes (particle content) files had to be given. Moreover because we use a non-linear 
gauge fixing condition|21j, the FeynArts default generic file had to be extended. 
This brings us to the issue of the gauge- fixing. We use a generalised non-linear gauge pUj 
adapted to the minimal supersymmetric model. The gauge fixing writes 

Cgf = -^lid^-iea-f^-igcwPZ^)^-^ + ^wkv + Sh + i:;H + ikx3)x''\^ 



-{d.Z + ^ZT^iv + ih + ^H)x,Y - (9.7) . 
2Ciy 

h and H are the CP-even physical Higgses, with h denoting the lightest. 7 is the photon 
field and the masses of the charged and neutral weak bosons are related through Mw = 
MzCw- The x's are the Goldstone fields. The non-linear gauge fixing parameters are a, 
(3, 6, u, k, e and 7. The ^ are the usual Feynman parameters. In our implementation the 
latter are set to ^ = 1 not only to avoid very large expressions due to the "longitudinal" 
modes of the gauge bosons but most importantly so that high rank tensors for the loop 
integrals are not needed. Gauge parameter independence which is a non trivial check on 
the result of the calculation can be made through the non-linear gauge fixing terms. In 
many instances a particular choice of the non-linear gauge parameter may prove much 
more judicious than another. For the case at hand, a = 1, preserves t/(l)em gauge 
invariance which explains the vanishing of the W~^x~'l vertex. This will prove crucial for 
the calculation of XiXi ~^ ^7 ■ 

This brinprs 11s tn thp imnlpmpntfltinn of thp loon intppTflls and thpir iisp in thp most prpripral 



application to radiative corrections in SUSY both for the colhders, for indirect detection 
and improvement of the rehc density calculation beyond tree-level. In LoopToolsj21] for 
example, the tensor loop integrals are reduced recursively to a set of scalar integrals by 
essentially following the Passarino-Veltman procedure [3 IJ. The reduction involves solving 
a set of equations that brings in the inverse Gram determinant^. Although for applications 
to the colliders the latter only vanishes for exceptional points in phase space, for the 
indirect detection calculation of tensor integrals involving annihilating LSP's with small 
relative velocity f , the Gram determinant is of order f ^. Therefore it vanishes exactly for 
V = or can get extremely small slightly away from this value rendering the calculation 
highly instable. In the Appendix we show how we dealt with this problem in an automatic 
implementation. In a nut-shell, we have used a segmentation procedure based on the fact 
that when some momenta are dependent like what occurs with f = 0, a A^-point function 
writes as a sum of A^ — 1 point functions. This also applies to the tensorial structures. 
This observation is not new (see for example E]) and has been used mostly in hand 
calculations. Some aspects of it may be remotely related to [S3]. The scheme also allows 
an expansion away from exactly vanishing Gram determinants. 




d) e) f) 

Figure 1: Typical classes of diagrams common to both XiXi ~^ 11 '^iT'd XiXi ~^ ^1 ■ For 
the former, the same "flavour" circulates in the loops and we do not have any mixing as 
in a) and d). Diagrams with Goldstone bosons are not shown. In the heavy wino and 
higgsino limit, a) is the dominant diagram. Diagrams d) and f) with H,h exchange do 
not contribute for v = 0. 

A selection of diagrams that contribute to both XiXi ~^ 11 XiXi ~^ ^1 is shown 
in Fig. n (see also E])- Diagrams of type a) in Fig. ^ are particularly important in 
the large wino and higgsino limit. In this limit the LSP and the internal chargino are 

""^For a recent overview of the problem with the Gram determinant, see \,V2\ . We will however present, 
for the 2^2 processes, a simple solution. 



of almost equal mass. If the W mass can be neglected this leads to a threshold singularity. 




Figure 2: An additional class of diagrams describing the Xil transition that only 

appear in the case of XiXi ~^ ^1 ■ ^ representative of the blob is the virtual correction. 
The counterterm contribution is shown also. 

Moving from XiXi ~^ 77 to XiXi ~^ -^7 brings mixing effects in the loops. Most diagrams 
can be derived from XiXi ~^ 77 • There is however an important class that is only 
present for XiXi ~^ ^1 ^is shown in Fig 121 This class of diagrams is missing in^J. It 
corresponds to the insertion of the Xi* ~^ Xi7 transition^. In a general gauge, the virtual 
transition would be gauge dependent and not ultraviolet finite. To remedy both these 
problems requires the Xi7 counterterm which is generated starting from the (tree-level) 
XiXi^ vertex through a Z — 7 one-loop transition. It is well known that the latter is 
gauge dependent, see for example ^21^. The counterterm requires the field normalisation 
5zV.^fn\. This field renormalisation constant in fact also induces, like in the standard 
model, {H, h)Z'y vertices not present at tree- level. This induced vertices are also needed 
for the class of diagrams shown in Fig. ^ in particular those with (H,h) exchange. The 
full set of counterterms for XiXi ^1 is in fact obtained from the tree-level XiXi ^ 
replacing a Z by a photon and inserting the ^Z^^ renormalisation constant. However, it 
is known that this renormalisation constant vanishes for a = We have checked this 

property explicitly with our code. After this check has been made, a = 1 was set, since 
it considerably reduces the number of diagrams and most importantly allows to discard 
all the counterterm contributions. Further gauge parameter independence of the result 
for XiXi ~^ ^1 "was checked by varying the other non-linear gauge parameters that enter 
the calculation, namely (3, 6 and a). When discussing our results for XiXi ^7 W6 will 
weight the effect of the new class of diagrams shown in Fig. |21 against that of the full 
contribution, in doing so we will specialise to a = 1. As we will see these diagrams give 
a non negligible contribution especially for the Higgsino case. 

The application to XiXi ~^ 99 is rather straightforward. This process confirms that our 
code handles the colour summation correctly. Keeping only one flavour of quark with 
charge Qf, XiXi ~^ 99 can be derived from XiXi ~^ 77 through {NcQfaY 2a'j. {Nc = 3 
is the number of colours, a is the electromagnetic fine structure constant and as the QCD 
equivalent). 

^The radiative neutralino decay x'j Xil is calculated inpi]. 



3 Results and comparisons 



We first check that our results are ultraviolet finite by changing the numerical value of 
the parameter 1/e that controls a possible ultraviolet divergence, e = 4 — n, where n is 
the dimensionality of space. We also check for gauge parameter dependence by varying 
the non-linear gauge parameters, namely <5, 6 and uj for XiXi ~^ 77 and /?, 5 and Cj with 
(5 = 1 for XiXi ~* ^7 • These checks are carried in double precision and show that, for 
all points we have studied, the results are consistent up to 13 digits. It is important to 
maintain the relation M^y = cy/Mz- If these parameters are taken as independent the 
gauge parameter independence is lost. We first discuss our results for 6 representative 
scenarios that we think are good checks on different parts of the calculations and also 
because they reveal the most important characteristics of these cross sections. Moreover 
these scenarios also serve to perform tuned comparisons against codes that are publicly 
available. To achieve this it is best to feed the codes the same parameters. Comparisons 
that use as input high-scale values for some SUSY parameters that are run down through 
some Renormalisation Group Equation, RGE, package often need to specify an interface. 
Moreover most often the RGE codes are updated and one does not always have access 
to the same version to perform a tuned comparison. For all these scenarios the input pa- 
rameters are defined at the electroweak scale and are: Mi the U{1) gaugino mass, M2 the 
SU{2) counterpart, the Higgsino "mass". Ma the pseudoscalar mass and the com- 
mon sfermion mass, tan/? is set to 10. The sfermion trilinear parameter Af is set to zero 
for all sfermions but the stop, depending on the mass of the latter. Our Higgs masses here 
are tree-level Higgs masses, so we avoided points too close to any Higgs resonance and the 
issue of the implementation of the width. When our code for XiXi ~^ 11 ~^ ^1 ^^id 

XiXi ~^ 99 will be fully incorporated within micrOMEGAs, corrected Higgs masses and 
mixing angles will be properly implemented in a gauge invariant manner following |35j. 
This improved implementation is of relevance only for f 7^ since the CP-even Higgses 
do not contribute when the neutralinos are at rest. When we refer to cross sections this 
would in fact refer to the cross section times the relative velocity v, av expressed in cm^/s. 
In terms of v, the total invariant mass of the neutralinos is s = 4m?o/(l — 
The six scenarios have been chosen so as to represent different properties, like the gaug- 
ino/higgsino content for different masses of the neutralino. We made no attempt whatso- 
ever to pick up points that lead to a good relic abundance in accord with WMAPPP. This 
said, for each point, we give the corresponding values of the relic density, extracted from 
micrOMEGAs. One should however keep in mind that different unconventional histories of 
the Universe^ could alter the usual thermal prediction. 

• Scenario 1: "Sugra". This reproduces a typical output from so-called mSUGRA 
scenarios, although the latter would not produce a common sfermion mass. The 
neutralino is mostly bino with mass around 200GeV. The lightest chargino is a 
wino. 

• Scenario 2: "nSugra". The neutralino is quite light, about lOOGeV and it is essen- 
tially bino. Here the mSugra relation does not hold, rather M2 = 4Mi. 

■^A few possibilities are described in|3ti|. 



• Scenario 3: "higgsino 1" . The neutralino is a light higgsino of about 200GeV. The 
hghtest chargino has a mass about 6GeV away. 

• Scenario 4: "higgsino 2" . The neutrahno is a heavy higgsino of about 4TeV. It is 
quite degenerate with the hghtest higgsino-hke chargino. The mass diference is 
about O.lGeV. 

• Scenario 5: "wino 1" . The neutrahno and hghtest chargino are hght, about 200GeV. 
The mass difference on the other hand is extremely small O.OlGeV. 

• Scenario 6: "wino 2" . This is like the previous example but for TeV masses. The 
LSP is a wino of mass 4TeV completely degenerate with the chargino. 

Table shows that the nature of the LSP and its mass critically determine its self- 
annihilation cross section to 77 and Z7. The results for the different scenarios vary 
by 6 orders of magnitude, especially for v = 0. The bino-like LSP gives far too small 
cross sections that are unlikely to be observed as a 7-ray line. The largest cross sections 
XiXi ~^ 77 cind XiXi ~^ ^7 fo^^ ^ LSP mass up to 4TeV are found for the wino-like LSP. 
Moreover in this case the signal is almost an order of magnitude stronger in Z7 than in 
77, however the two lines even for M^o ~ 200GeV are only lOGeV away, even before 
any smearing is taken into account. For the wino case the contribution of the — >■ x57 
transition to XiXi ~^ -2^7 is negligible for the two scenarios we display in Table [U This is 
not the case of the higgsino scenarios (nor the bino-like for that matter) where this con- 
tribution could amount to a correction of more than 30%. It is also interesting to note, 
see later, that for the very heavy wino scenario the cross section drops very quickly as we 
increase the velocity. We will study the wino and higgsino case in more detail below. 
For XiXi ~^ 99 ? the LSP composition does not show dramatic differences in the cross 
sections. The largest are however found for a light wino and a light higgsino of 200GeV. 

Let us now turn to the comparisons, essentially for v = 0, with the codes PLATON and 
DarkSUSY. For this fine-tuned comparison we have taken the same masses for the quarks 
and Mz as in DarkSUSY ^2] as well as for the electromagnetic and strong coupling. On 
the other hand we imposed Mw = MzCw- Taking for example, Mw = 80.33GeV, with 
all other parameters as in Table [T] not only gives gauge parameter dependent results, but 
in the wino case the LSP would turn out to be the chargino. 

Table Q shows that our results (for f = 0) agree perfectly with those of PLATONdml as 
concerns XiXi ~^ 77 as well as with PLATONdmg as concerns XiXi ~^ 99 ■ PLATONgrel 
also perfectly confirms our results for XiXi ~^ 99 for v = 0.5. Excellent agreement with 
DarkSUSY is also observed (at f = 0) for XiXi ~^ 11 and XiXi ~^ 99 ■ To compare with 
the results of DarkSUSY for XiXi ~^ ^7 do not consider the contribution from the 
Xi ^ Xil "insertion" . With this restriction we find exactly the same results as DarkSUSY 
in the case of the wino and higgsino but not in the case of the bino, where our results are 
about 30% higher. However as pointed out earlier, the cross sections in the bino case are 
tiny. 

The effect of the new contribution from Xi ~^ X1I foi' both v = and v = 0.5 is not 
noticeable when the neutralino is a pure wino, but it can be important in the higgsino 
case where the total contribution is also large. The new contribution brings in a relatively 
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Table 1: Results of our calculation both at v — and v — 0.5 for XiXi ~^ 77 > XiXi ~^ 
and XiXi -^7 O'lT'd comparison with the codes o/PLATON anc? DarkSUSY. Forx^Xi ~^ 
, 'full" refers to the compete set of diagrams. "Part" refers to excluding the Xi'Xi7 

insertion. Inputs are at the electroweak scale and are expressed in TeV. tan/? = 10. 

Af = apart from At = —300GeV for m? = 0.8 TeV and At = for m? = 10 TeV. 

We have taken the DarkSUSY inputs, with Mz = 91.187 GeV and = 0.2319 (but Mw = 

MzCw)- The quark masses are rut — 175 GeV, rub —5GeV,mu — 56MeV, rud — 99MeV, 

rus = 199Mey, rUc = 1.35GeV. = 127.942 and as = 0.117. The relic abundance Qh'^, 

extracted from micrOMEGAs, is also given for completeness. 



small correction in the bino case, where the cross sections are tiny anyhow. 
Results for XiXi ~^ 77 XiXi ~^ ^7 for = 0.5 have been computed here for the first 
time and can be relevant for the computation of the relic density in some regions of the 
parameter space. 



4 The wino and higgsino limits 




Figure 3: Dependence of the 77 and Z7 cross sections as a function of the LSP neutralino 
mass M^o. For the higgsino case we take M2 = 2Mi = 4 lO^TeV and for the wino 
;x = Ml = 2 10^ TeV. In both cases tan (3 = 10, Ma = 100 TeV and all sfermions 4 10^ TeV. 
We use Mw = MzCw, with sw = 0.473, Mz = 91.1884^61/, = 127.9. 

The results of Tabled make it clear that most interesting scenarios for the monochro- 
matic 7 ray line signals are of a wino and higgsino type even when the LSP has a mass 
of about 2Mz. Fig. 01 shows the dependence of the XiXi ~^ 77 and XiXi ~^ ^7 cross 
sections at f = as a function of the LSP mass in the case of a wino and a higgsino LSP. 
The mass of the LSP is in the range 70GeV to lOOTeV. In fact masses below lOOGeV 
may be excluded by LEP2 but it is interesting to see how the cross sections grow past the 
lOOGeV mass to stabilise around a plateau. The masses of the other supersymmetric par- 
ticles are taken extremely heavy here. Note that in the higgsino case XiXi ~^ ^1 shows 
much more structure. The peak cross section is much more pronounced before the cross 
section decreases and reaches a plateau of the same order as XiXi ~^ 77 • The wino cross 
section for XiXi ~^ ^1 is the largest of all and is almost an order of magnitude larger 
than XiXi ~^ 11 ^i-nd two orders of magnitude larger than in the higgsino case. 

It is also interesting to see how the plateau is reached for a fixed mass of the wino and 
hip'P'sinn T;SP. or rather a fixed value of Mn and denendiupr on the comnosition of the 




Z 7 (wino) 



7 y (wino) 



10 



y y (higgsino) 



10^ 



10 



10 

li(TeV 



01 


io" 






E 
o 


10" 


> 
b 


10" 




io" 




io" 




io" 



-26 



28 



29 



31 







Z 7 (wino) 








y y (wino) 










y y (higgsino) 




/ - - - 

* 









1 1 









10 



10" 



10 



li(TeV 



10" 



Figure 4: The first figure shows the dependence on Mi for fi{M2) fixed in the higgsino 
(wino) limit for 77 and Z'-f. The values for the SUSY parameters are: tan/5 = 10, 
= 4 lO^GeV, A = 0, Ma = 100 TeK /i = bOTeV, M2 = 2Mi m the Mggsmo case. 
M2 = 10 TeV and fi = Mi in the wino case. In the figure at the bottom, the only parameters 
that are changed are fi = 10 TeV in the higgsino case and tan/? = 2, M2 = 500 GeV in the 
wino case. The SM parameters are as in Fig. O 



LSP. We therefore fix these two values and vary the other supersymmetric parameters 
of the neutrahno sector. The behaviour of the cross section as we vary these parameters 
is shown in Fig. EJ For the wino case one can see that once Mi, fi > M2, and therefore 
the LSP is mostly wino, the asymptotic values are reached abruptly especially for the 
case of a wino of lOTeV and large tan f5. Below this transition, the cross sections have a 
smooth behaviour. In the higgsino limit, a fast transition occurs once Mi,M2 > fi but 
past this threshold there is still a smooth and slow increase of the cross sections before 
the asymptotic values are reached. 

Most of this behaviour can, in fact, be recovered through simple analytical expressions 
that serve also as a further check on our results and the accuracy of the calculation in these 
extreme scenarios. It had been observed ITT] that when the LSP is heavy, much heavier 
that the 14^-boson, the cross sections (times velocity) for XiXi ~^ 77 and XiXi ~^ -2^7 tend 
to an asymptotic value that can be computed from the dominant contribution, that of 
diagram a) of Fig. ^ The limiting behaviour can be easily understood from the fact that 
in the heavy mass limit, the annihilating LSP neutralino and the chargino are degenerate 
with a mass much larger than the weak boson mass. This develops a threshold singularity 
like what one finds in QED, although here acts as a regulator. Another important 
factor that measures how the asymptotic values are reached is the deviation from exact 
degeneracy between the LSP neutralino and the lightest chargino given by their mass 
difference, 6m, 6m = m-+ — m^ojH]. This scenario has been revisited in a series of 
excellent papers^!] where it has been shown how the one- loop calculation in these cases 
need to be improved through a non-relativistic treatment. Our aim here, in the rest of this 
section, is to see how our one-loop results can be made to match with the non-perturbative 
treatment. This paves the way to an implementation in a code for indirect detection that 
can be used in all generality like what we have started to do in micrOMEGAs. 
First, we will show how our one-loop results effectively capture the behaviour of the cross 
sections in these scenarios and how the asymptotic value in the case of a wino is reached 
dramatically fast. 

In the higgsino limit, /j, -C Mi, M2. We will also take M2 = 2Mi = 2Ms and consider also 
the large tan/5 (in fact tan/? > 2 suffices), /i will be taken positive. 
In the wino limit, M2 <C fJ-jMi. We will also take fi = Mi = Ms and large tan/?. The 
(tree-level) mass difference in the higgsino, 6m^, and wino limit, 6m'^, write 



We see that in the wino case the mass difference scales like l/M|jT^^ compared to the 
I/M5 in the higgsino case. In these configurations the cross sections are well approximated 
by a^'^'^ V in the higgsino case and a^'^'"' v in the wino case {V = Z, 7) which are the 
results of the dominant diagrams of Fig. Q-a), 



''it is important to note that it is essential to have AIw = MzCw, otherwise we could get a mass 
difference oc 1/Ma. 



al^'^v I 1 + 



2M|M2 
^ Ml tan (3^ 



(2) 



where the asymptotic values {5m = 0) are given by 



4 

an 



16(t2^'~^ V ~ 1.6 IQ-^^cmys, 



[1/2-8 



2 \2 

w) 



„2 2 



a. 



■2^'^ V ~ 0.8 IQ-^^cm^s, 



c 



Q -W 77,10 
^ „2 '^oo 



V ~ lO^^^cmVs 



(3) 
(4) 

(5) 



We have verified that our code including the complete contributions agrees extremely 
well with these approximation for the cross sections, Eq. |21 and that moreover in the wino 
case the asymptotic values, Eq. 01 are reached very fast due to very small 6m. This is 
also exemplified in Fig. HJ 

As demonstrated in ^3] the one-loop treatment of the threshold singularity that is 
responsible for the behaviour of these cross sections in the higgsino and wino regime at 
high LSP mass is not adequate an breaks unitarity. The non-perturbative non-relativistic 
approach of Ref. ||14j not only improves on the calculation but it also unravels the forma- 
tion of bound states that drastically enhance the annihilation cross sections for specific 
combinations of masses. Fig. shows the effect of such resonances and the departure of 
the cross section XiXi ~^ 11 from a full one-loop treatment as the mass of the higgsino 
LSP increases. The "resonance" curve is based on the use of a fitting function as given 
in|14j ^. All other particles are taken heavy apart from the higgsino mass parameter /i and 
Mi = Ms = M2/2 = 25TeV. For the whole range of M^o ~ /i we have 5m = O.lGeV. The 
figure also shows the value of the approximate one-loop result as given for the higgsino 
limit in Eq. El together with the full one-loop treatment based on our calculation. The 
resonance formation, here around 6TeV, brings an enhancement factor of more than 4 
orders of magnitude. On the other hand departure from the full one-loop calculation is 
of relevance only for M^o masses around ITeV. The insert in Fig. (HI shows in more detail 
the comparison of the full calculation compared to the approximate result for the smaller 
higgsino LSP masses, well before the resonance effects settle in. For M^o > 600GeV the 
approximation is very good, only around M^o ~ 200GeV, the full calculation captures the 
effect of other contributions like those of Fig. ^,e). For this particular case it looks like 
a good matching between the full one-loop result and the non perturbative one should 

^We thank J. Hisano and M. Nojiri for confirming that Eq. 61 of hep- ph/041 2403, should be squared 
and that the entry i = 0, 7 = 1 in Table 1 of that paper is 10 times smaller. 




Figure 5: Comparison, for XiXi — * 77 os a function M^o, between our prediction for the 
full one-loop calculation (continuous line, black), analytical one-loop expressions based on 
the approximation of Eq. (dotted, blue) and the non-perturbative prediction based on 
the fitting functions including resonances'^^ (dashed, red). The input susy parameters are 
M2 = 2Mi = 50TeV, tan/3 = 10, A = 0, mj = Ma = lOOTeV. The SM parameters are 
as in Fig. The insert is a close-up for small M^o . 



be made at around M^o = 400GeV. A possible strategy for choosing this matching point 
would have to rely on the knowledge of both the full one- loop result, the approximate one- 
loop result as given Eq.|21and the non-perturbative results based on the fitting functions of 
Ref . |I4j . This would, of course, only be carried out in the limit of almost pure higgsino or 
wino. We would then have to compare the three results. To revert to the non-perturbative 
regime means that the full one-loop and the approximate one-loop agree fairly well and 
are quite different from the non-perturbative regime. If, on the other hand, these two 
one-loop results differ sensibly this means that one is not quite in the asymptotic region 
and that we might be missing some one-loop contributions. If this is the case one should 
also expect the higher order effects computed for the threshold region to be small so that 
the non-perturbative result and the approximate one-loop are very similar. Of course, 
as shown in the example of Fig. these differences in the higgsino region, compared to 
taking the perturbative parameterisation, are rather small compared to the uncertainty 
that is inherent in the astrophysics part of the prediction of the gamma ray line. For the 
wino, as we saw, the transition to the asymptotic value is rather drastic especially for 
TeV LSP's, therefore one should quickly capture the non-perturbative regime. Especially 
in this case one should also revert to a one-loop use of the chargino-LSP mass difference. 



5 Conclusions 



There has been a flurry of activity in the last few years in the search of dark matter 
with, among other strategies, many experiments dedicated to the indirect searches of 
dark matter in particular the gamma ray signal. The mono-energetic gamma ray line 
signal constitutes a clean signature. The improvement in coverage and accuracy of the 
measurements should be matched by precise theoretical calculations that should be pub- 
licly available through general purpose codes. In this paper we have provided a new 
calculation for XiXi ~^ 77 ^ind XiXi ~^ ^1 fo^' the annihilation of the supersymmetric 
dark matter candidate and rederived as a bonus the XiXi ~^ 99 rate. These calculations 
have been made both for small (zero) relative velocity of the neutralinos as adequate for 
annihilation in the halo of our galaxy for example, but also for velocities that would be 
needed for the contribution of these channels in a precise derivation of the relic density. 
For XiXi 77 and XiXi ~^ ^1 at f 7^ as would be needed for an improved relic density 
prediction, these results are new. For XiXi ~^ 77 three jHl 1101 full one-loop calculations 
performed for v = have already been performed. We have performed a tuned compar- 
ison with the results of DarkSUSY^2] and PLATON^n] and have found perfect agreement. 
The calculation of XiXi ^7 is trickier and can not just be deduced from XiXi ~^ 77 • 
Until now there has been only one calculation^!] of this process. The latter has missed 
some contributions that may not always be negligible. Comparisons of our results with 
the previous ones without these contributions are quite good for scenarios where the cross 
section is not small. In this paper we have not made an attempt to fold in with the 
astrophysics part that involves, for example, the halo proflle but concentrated on the 
particle physics part which must be an unambiguous prediction. To pave the way for an 
implementation in micrOMEGAs[15^ we felt it was important to critically review the large 
mass higgsino and wino LSP scenarios especially that the latter give large cross sections. 
As shown in ^3] one needs to go beyond the one-loop treatment in this regime. We have 
argued how one could match the full one- loop treatment with the non-perturbative result. 

Another important aspect of this paper is the way all these calculations have been per- 
formed in a unifled manner and the techniques that we used. These processes are the flrst 
application of a code for the calculation of one-loop processes in supersymmetry both at 
the colliders and for astrophysics/relic density calculations that require also a new way of 
dealing with the reduction of the tensor integrals. The calculations are performed with 
the help of an automatised code that allows gauge parameter dependence checks to be 
performed. The use of a generalised non-linear gauge is crucial. 
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A Appendix: Segmentation of loop integrals 

The tensor integral of rank M corresponding to a A^-point graph, {M, A^}, that we en- 
counter in the general non-linear gauge but with Feynman parameters ^ = 1 are such 
that M < N. For the evaluation of 2 ^ 2 processes A^ = 4 is the maximum value and 
corresponds to the box. The general tensor integral writes as 

M 

where 



Di = {l + SiY-M^, s, = Y.p,, so = 0. (A. 2) 

Mi are the internal masses, pi the incoming momenta and / the loop momentum. 

The A^-point scalar integrals correspond to M = 0. All higher rank tensors for a 
A^-point function, M > 1, can be deduced recursively from the knowledge of the A^-point 
(and lower) scalar integrals. In Looptools|24j this is based on the Passarino-Veltman 
algorithm |31j. In Grace-loop the implementation is outlined in Ref. [21]. The tensor 
reduction involves solving, recursively, a system of equations which explicitly requires 
the evaluation of the Gram determinant: DetG{pi,p2,P3) = DetGi23 = DetpiPj. For 
special kinematics the latter vanishes or can get very small, leading to severe numerical 
instability. This special kinematics for the general 2 — >^ 2 process one encounters in high 
energy occurs for exceptional points in phase space, for instance in extremely forward 
regions and most generally the weight of this contribution may be dismissed. For the case 
at hand, when the two neutralinos are at rest, or with extremely low relative velocity, 
the Gram determinant vanishes for all points because the incoming neutralinos have the 
same momentum and can not be considered independent. This is exactly what occurs in 
our case. Indeed here, the box diagrams for XiXi ~^ ^1 have the Gram determinant 



^etG = M|o.^^^— ^(1-.^), -^ = 4^(1 --74), (A. 3) 



w is the relative velocity and 9 is the scattering angle. For XiXi — > 77 , = 0. In our 
application the sub-determinant with the incoming LSP neutralinos is responsible for the 
vanishing of the box Gram determinant for all angles: 



This also means that the reduction of the tensor integrals for triangles with the two LSP 
as external legs, needs special treatment. Such triangles are of the type Fig. l.e) (but not 
Fig.l d)). 

We therefore need to implement a routine for cases when the determinant vanishes 
due to the fact that the momenta Si,i ^ are not independent. There are a few ways 
of dealing with the tensor integrals when the Gram determinant is exactly zero [221 ISS] • 
Sometimes, the problem is even avoided by the grouping of terms so that this spurious 
inverse determinant cancels out. Our aim was to find, at least for 2^2 process, an 
efficient method that can, not only be easily automated, but also calls most of the existing 
routines that are present for a general purpose algorithm designed for non zero Gram 
determinants. Take the box for example. Observe that (in any n dimension), in most 
generality, we may write for any given pair of constants a, (3 



a - (3 + (a + ^ - 1)- 



1 



A + 2l.{s^- asi- (3s2) 
A = {sl-M',)~a{sl~M!)-(3{sl-Ml)~{a + (3-l)Ml (A. 5) 



Obviously if S3 = asi + /3s2 and hence the momenta are linearly dependent, the box 
splits into a sum of triangles. We will refer to this as segmentation. This segmentation 
is independent of the tensor structure. This means that the reduction of the tensor box 
with zero Gram determinant amounts to a tensor reduction for triangle with a non-zero 
Gram determinant for which one uses the usual procedure and hence uses the general 
library. Observe that if a = or /5 = 0, there are three segments instead of four. The 
missing segment, triangle integral, does in fact have a zero Gram determinant. Therefore 
when one approaches the zero of the Gram determinant of the box in these specific cases, 
a ~ for example, a will be numerically very small but non zero. A numerical instability 
could still develop due now to the Gram determinant of the associated triangle. These 
"algebraic" zeros could be missed at the numerical level. This can again lead to (less 
severe) numerical instabilities due to the reduction of the associated tensor. In this case, 
these triangles are segmented even further into two-point functions, following the same 
recipe. This wav their contribution is negligible even at the numerical, automatic, level. 



In any case, since we also encounter (tensor) triangle diagrams (Fig. l.e)) that have a 
vanishing Gram determinant, we have also included a segmentation that also applies to 
the tensor triangles using the same trick as for the boxes. 

There is an important observation to make. The segmentation of a tensor of rank M for 
the iV-point function, {M, N}, amounts to applying the tensor reduction on {M, — 1}. 
If M = iV, after segmentation one would need a library for {N,N — 1}. These are 
not supplied by default in the general libraries. These libraries would then need to be 
extended. Reduction of — 1-point function tensor integrals are much more compact 
and easier than for A^-point functions. This said for the case at hand, and for that matter 
any relic density calculation where the LSP is a neutralino, these highest rank tensors are 
not needed. It is easy to show that the highest rank tensors for 2 — > 2 only occurs when 
all the external particles are bosons. In our case, for the box, one has Mmax = 3. For 

The choice of the momenta circulating in the loop, 81,82,83, is not unique and depends 
on the particular graph. If DetG(si, S2, S3) = we first form all three sub-determinants 
DetG{si, Sj) and take the couple Si,Sj that corresponds to Max \Det{si, Sj)\. Then the 
third Sk is distributed in the basis Sj, sj and the corresponding a and (3 are read. In fact, 
suppose DetG{8i, 82) 7^ 0, then it is revealing to always write 



S3 = asi + f3s2 + £t with 
si.St = S2.eT = meaning 
£t,h = e^apssis^t^. (A. 6) 

Et is a vector that is orthogonal to both si and S2 that is easily reconstructed once a and 
P are 



sjss.Si - 81.8282.83 o ( N ( \ ^\ 

^ = nTTv ^ ' 13 = oi[si^ S2). A. 7 

DetG{8i, 82) 



This construction makes it clear that 



DetG{8i, 82, 83) = elDetG{8i, 82). (A. 8) 

This shows, in a most transparent manner, that the determinant vanishes when — 0. 
This can occur when the components of this vector vanish, et^ 0, and therefore S3 is not 
an independent vector as the case of this paper for v = 0. It also occurs, a point which 
is often overlooked, when Et is light-like. However the orthogonality constraint means 
that the other vectors si,S2 are space- like. Therefore this case does not occur for real 
particles that are time-hke, and hence does not occur for our 2 — > 2 process. It will be 
shown, in a separate publication, that when Et is light-like a segmentation is still possible. 
The algorithm can also be improved by expanding around DetG{si, S2, S3) = 0. We will 
come back to the details of this issue in a future publication. Note that for v = 0.5 and 
for all three processes studied in this paper the standard reduction formalism, without 

seprnentatinn . was iiserl. 
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